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ABSOLUTE PROPERTIES OF THE SPOTTED ECLIPSING BINARY STAR CV BOOTIS 

GUILLERMO T0RRES\ LUIZ PAULO R. VAZ^, AND ClAUD H. SANDBERG LACY'' 
Draft version September 12, 2008 

ABSTRACT 

We present new V^-band differential brightness measurements as well as new radial-velocity mea- 
surements of the detached, circular, 0.84-day period, double-lined eclipsing binary system CV Boo. 
These data along with other observations from the literature are combined to derive improved absolute 
dimensions of the stars for the purpose of testing various aspects of theoretical modeling. Despite 
complications from intrinsic variability we detect in the system, and despite the rapid rotation of the 
components, we are able to determine the absolute masses and radii to better than 1.3% and 2%, 
respectively. We obtain Ma = 1.032 ± 0.013Mq and i?B = 1.262 ± 0.023 Rq for the hotter, larger, 
and more massive primary (star A), and Mb — 0.968 ± O.O12M0 and R-q — 1.173 ± 0.023 R© for the 
secondary. The estimated effective temperatures are 5760±150K and 5670±150K. The intrinsic vari- 
ability with a period ^1% shorter than the orbital period is interpreted as being due to modulation by 
spots on one or both components. This implies that the spotted star(s) must be rotating faster than 
the synchronous rate, which disagrees with predictions from current tidal evolution models according 
to which both stars should be synchronized. We also find that the radius of the secondary is larger 
than expected from stellar evolution calculations by ~10%, a discrepancy also seen in other (mostly 
lower- mass and active) eclipsing binaries. We estimate the age of the system to be approximately 
9 Gyr. Both components are near the end of their main-sequence phase, and the primary may have 
started the shell hydrogen-burning stage. 

Subject headings: binaries: eclipsing — stars: evolution — stars: fundamental parameters — stars: 
individual (CV Boo) — stars: spots 



1. INTRODUCTION 

> ■ CV Boo (= BD +37 2641 = GSC 2570 0843; a = 
' 15^ 26"^ 19554, S = 4-36° 58' 53'.'4, J2000.0; V « 10.8, 
OO . SpT = G3V) was discovered as a p ossible e clipsi ng bi- 
CN ■ nary star bv iPeniche et al.l (|1985f ). iBuschI ()1985f ) con- 
firmed it as an eclipsing binary of type EA and found its 
Q\ . period to be 0.8469935 days. In his last p ublished paper 
' a study of 4 lower main sequence binaries. iPopped (|200Gl ) 
OO , determined a spectroscopic orbit for CV Boo. Popper 
. was pessimistic about the prospects for determining ac- 
L| ■ curate absolute properties of them because "It appears 
. ^ unlikely that definitive photometry will be obtained for 
these stars, partly because of intrinsic variability." Re- 
$H cently, a light curv e and radial vel ocity study of the sys- 
5^ [ tem were done by iNelsonI (j2004b| ) , resulting in the first 
estimates of its absolute properties. 

The parameters of CV Boo make it potentially in- 
teresting as the most evolved system among the well- 
studied double-lined eclipsing binaries with components 
near 1 Mq (see Figure[T]), a regime where some discrepan- 
cies with theoretical models have been pointed out. We 
describe in the following our extensive new photometric 
and spectroscopic observations of the object intended to 
improve our knowledge of the system. The presence of 
starspots does in fact limit somewhat our ability to deter- 
mine highly accurate absolute properties for this binary 
star, but the results are still accurate enough for mean- 
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Fig. 1. — Main-sequence eclipsing binaries in mass range of 
CV Boo with accurate determinations of their absolute properties 
(masses a n d rad ii good to better than ~2%). Data are taken from 
[AndersenI 1119911 ) and updates from the literature. Primary and 
secondary components are connected with solid lines. CV Boo is 
represented with open circles. The dashed line shows the solar- 
metallicity zero-age main sequence from the models by lYi et al.l 
II2001I '). for reference. 



ingful tests of current stellar models. As we describe 
here, CV Boo contributes significantly to the body of ev- 
idence concerning the differences with theory mentioned 
above. 
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2. OBSERVATIONS AND REDUCTIONS 
2.1. Differential and absolute photometry 

New differential brightness measurements of CV Boo 
were obtained with the faciUties available at the Kimpel 
Observatory (ursa.uark.edu). They consist of a Meade 
10-inch f/6.3 LX-200 telescope with a Santa Barbara In- 
struments Group ST8 CCD camera (binned 2 x 2 to pro- 
duce 765x510 pixel images with 2.3 arcsec square pixels) 
inside a Technical Innovations Robo-Dome, controlled 
automatically by an Apple Macintosh G4 computer. The 
observatory is located on top of Kimpel Hall on the 
Fayetteville campus of the University of Arkansas, with 
the control room directly beneath the observatory inside 
the building. Sixty-second exposures through a Bessell V 
filter (2.0 mm of GG 495 and 3.0 mm of BG 39) were read 
out and downloaded by ImageGrabber (camera control 
software written by J. Sabby) to the control computer 
over a 30-second interval, and then the next exposure 
was begun. The observing cadence was therefore about 
90s per observation. The variable star would frequently 
be monitored continuously for 4-8 hours. CV Boo was 
observed on 89 nights during parts of two observing sea- 
sons from 2001 December 1 to 2003 June 9. 

The images were analyzed by a virtual measuring en- 
gine application written by Lacy that flat-fielded the im- 
ages, automatically located the variable and comparison 
stars in the image, measured their brightnesses, sub- 
tracted the corresponding sky brightness, and corrected 
for the differences in airmass between the stars. Ex- 
tinction coefficients were determined nightly from the 
comparison star measurements. They averaged 0.20 
mag/airmass. CV Boo is also known as GSC 2570 0843. 
The comparison stars were GSC 2570 0511 ("comp", 
V = 10.26, as listed in the Tycho Catalogue), and 
GSC 2570 0423 ( "ck" ) . Both comparison stars are within 
8 arcmin of the variable star ("var"). The compari- 
son star magnitude differences (comp— ck) were constant 
at the level of 0.013 mag (standard deviation within 
a night), and 0.007 mag for the standard deviation of 
the nightly mean magnitude difference. The differen- 
tial magnitude (var— comp) of the variable star was ref- 
erenced only to the magnitude of the comparison star, 
comp. The resulting 6500 V-band magnitude differences 
(var— comp) are listed in Table[T]and plotted in Figure^ 
The typical precision of the variable star differential mag- 
nitudes is about 0.013 mag per measurement. 

In addition to our own, differential photometry of 
CV Boo was obtained by INelsonI (|2'004bf ) in V and 
Cousins / between 2003 March and June (253 and 265 
measurements, respectively). The comparison star was 
GSC 2570 0869, and the check star was GSC 2570 0511, 
which is the same star we used as the comparison. These 
observations are incorporated into our analysis below. 

Absolute photometry of CV Boo is available in the 
literature from several sources, and color indices can 
be used to estimate a mean effective temperature for 
the combined light of the system (assuming no inter- 
stellar reddening). The results are collected in Ta- 
ble [2l We used the color /temperature calibrations of 
iRamfrez fc MelendezI (^005) for dwarf stars for all but 
the Sloan g — r index; for tha t color we used the cali- 
bration of iGirardi et al.l (j2004l ). In all cases we assumed 
solar metallicity. The value of Johnson V is that listed 
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Fig. 2. — The top panel shows our y-band hght curve of CV Boo 
(+ symbols), consisting of 6500 points, together with the V (x) 
and / (•) light curves from Nelson (2004b), marked "AVn" and 
"A/n". Nelson's light curves are shifted as indicated, for clarity. 
Our theoretical solution without spots is overplotted (continuous 
grey lines; ^ 13.11 1. The lower panels show the O — C residuals from 
these fits, and in the upper left corner, the standard deviation for 
a single measurement. 

in the Tycho Catalogue with no uncertainty given there. 
We have assumed a conservative error of 0.10 mag for 
V. The temperature estimate from the Johnson B — V 
index uses the value of that index as listed in the Tycho 
Catalogue, with its listed error. The temperature values 
estimated in these ways agree quite well, except for the 
estimates from B—V and B^—Vt, which happen to have 
the largest formal errors. The weighted average of the 
7 estimates is 5706 ± 60 K, where the uncertainty does 
not account for possible systematic errors in the various 
calibrations. This color-index-based temperature is quite 
consistent with spectroscopic estimates discussed below, 
and this suggests that the interstellar reddening value, if 
any, is very small. 

CV Boo is ide ntified as a strong X-ray source in the 
ROSAT catalog (jVoges et al.||2000( ). This is presumably 
due to an active chromosphere/corona associated with 
its spot activity (see below). 

2.2. Ephemeris 

Photoelectric or CCD times of minimum light for 
CV Boo carried out over the past decade have been re- 
portcd by a number of authors fAgerer fc IIiibscheii l2002|. 
|2p03; Bakis ct al. 2003; Diethclm 2001; Dogr u et al.l 
20061; Hiibschcr 2005; Hiibsc her et all I200a l2006t 
Kim et al. .2006.; .Lacv..20q2.. .2003inMacieiewski fc Karskal 
200l INelsonI [2003. 120021 l2004al) . Additional times of 
eclipse including older visual and photographic measure- 
ments reaching back to 1957 were kindly provided by J. 
M. Kreiner (se elKrei ncr. Kim fc Nha 2000) or taken from 
the literature l|Locher 2005 ; Molik 2007 ). Separate least- 
squares fits to the 98 available primary and 50 secondary 
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minima yielded ephemerides with virtually the same pe- 
riod within the errors. A simultaneous fit to all minima 
was then performed assuming a circular orbit. Uncer- 
tainties were initially adopted as published, or assigned 
by iterations and by type of observation so as to achieve 
a reduced of unity. The resulting linear ephemeris is 

Min I (HJD) 2,452,321. 845322(50)+0.846993420(69)S. 

where the figures in parentheses represent the uncer- 
tainty in units of the last decimal place. No significant 
trends indicative of period changes are seen in the O—C 
residuals. A test solving for separate primary and sec- 
ondary epochs with a common period yielded a phase dif- 
ference between the eclipses of A0 = 0.49991 ± 0.00013. 
This is consistent with 0.5, supporting our earlier as- 
sumption of a circular orbit. 

2.3. Spectroscopy 

CV Boo was observed spectroscopically with an echelle 
instrument on the 1.5m Tillinghast reflector at the F. L. 
Whipple Observatory (Mt. Hopkins, Arizona). A total of 
66 spectra were gathered from 1991 June to 2005 April, 
each of which covers a single echelle order (45 A) cen- 
tered at 5188.5 A and was recorded using an intensified 
photon-counting Reticon detector. The strongest lines 
in this window are those of the Mg I h triplet. The re- 
solving power of this setup is A/AA « 35,000, and the 
observations have signal-to-noise ratios ranging from 13 
to 36 per resolution element of 8.5 kms~^. 

Radial velocities were obtained using the two- 
dimensional cross-cor relation algorithm TODCOR 
(jZucker fc Mazehl Il994l ). Templates for the cross cor- 
relations were selected from an extensive library of 
calculated spectra based on model atmospheres by R. L. 
Kurucz'* (see also iNordstrom et al.l Il994t iLatham et al.l 
I2OO2I ). These calculated spectra cover a wide range of 
effective temperatures (Toff), rotational velocities (usini 
when seen in projection), surface gravities (logg), and 
metallicities. Experience has shown that radial veloc- 
ities are largely insensitive to the surface gravity and 
metallicity adopted for the templates, as long as the tem- 
perature is chosen properly. Consequently, the optimum 
template for each star was determined from extensive 
grids of cross-correlations varying the temperature and 
the rotational velocity, seeking to maximize the average 
correlation weighted by the strength of each exposure. 
Solar metallicity was assumed. The results, interpolated 
to surface gravities of \ogg = 4.25 for both stars (see 
§[3]), are Tcff = 5800 K and wsinz — 73 kms~^ for the 
primary star, and Toff — 5650 K and wsini = 67 kms~^ 
for the secondary. Estimated uncertainties are 200 K 
and 10 kms~^ for the temperatures and projected 
rotational velocities, respectively. Template parameters 
near these values were selected for deriving the radial 
velocities. Typical uncertainties for the velocities are 
5.6 kms~^ for the primary and 5.9 kms~^ for the 
secondary, which are considerably worse than usual with 
this instrument because of the significant rotational 
broadening of both stars. 

The stability of the zero-point of our velocity system 
was monitored by means of exposures of the dusk and 

* Available at |http: //cf akuS . cf a. harvard ■ edu| 



dawn sky, and small run-to- run corrections were applied 
in the manner described bv iLathaml ()1992l ). Additional 
corrections for s ystematics were appli ed to the velocities 
as des cribed by ILatham et al.l ()1996f ) and [Torres et al.l 
()1997f ) to account for residual blending effects and the 
limited wavelength coverage of our spectra. These cor- 
rections are based on simulations with artificial compos- 
ite spectra processed with TODCOR in the same way 
as the real spectra. The final heliocentric velocities are 
fisted in Tabled 

The light ratio between the compone nts was estimated 
directly from the spectra following IZucker fc Mazehl 
()1994D . After corrections for systematics analogous to 
those described above, we obtain ^b/^a = 0.71 ± 0.04 
at the mean wavelength of our observations (5188.5 A), 
where we refer to star A as the more massive one (the pri- 
mary) and to the other as star B. This valu e is in reason- 
able agreement with estimates by 'Popped (|2000f ) based 
on the relative strength of the Na I D lines in CV Boo. 
Given that the stars have slightly different temperatures 
(see below), a small correction to the visual band was 
determined from synthetic spectra integrated over the V 
passband and the spectral window of our observations. 
The corrected value is (^b/^a)^ = 0.73 ± 0.04. 

Radial veloci ties for CV Bo o have been reported pre- 
viously also bv iPoppeil (|2000( ). who observed the star as 
part of his program focusing on binary systems of spec- 
tral type F to K. His 45 measurements from 1988 Febru- 
ary to 1997 June with the Hamilton spectrometer at the 
Lick Observatory partially overlap in time with ours, and 
are of excellent quality. However, they require a number 
of adjustments before they can be combined with ours. 
One of these adjustments has to do with corrections he 
applied to his raw velocities. The raw velocities (which 
he referred to as "Observed" ) were reported for CV Boo 
alongside "Orbital" velocities which differ from the raw 
ones by the application of two corrections. The first is 
analogous to the corrections for systematic effects we ap- 
plied to our own velocities, and was derived in a similar 
man ner using synthet ic binary spectra (for details see 
Popper fc Jeong|[l994f ). The second correction accounts 
for distortions and mu tual irradiation in the close orbit, 
and was computed by 'Popper' (| 200(3) using the formal- 
ism developed bylWilson (1990) as implemented in the 
Wilson-Devinney (WD) program that we also use below, 
and added to the velocities. Given our plan to use the 
WD program to combine the light curves with the veloc- 
ity measurements in a si multaneous solu tion, the latter 
corrections in the data by iPoppeH (|2000f ) need to be re- 
moved prior to use or they would be applied twice. We 
estimated these corrections from a preliminary solution 
with WD (presumably emulating Popper's procedure), 
and applied them with the opposite sign to the "Orbital" 
velocities. These corrections are no larger than 1 kms^^, 
which is smaller than the formal uncertainties in the ve- 
locities (2.7 kms~^ for the primary and 2.1 kms^^ for 
the secondary, from the residuals of preliminary fits). 

A second adjustment we found necessary to apply to 
Popper's measurements is to correct for an offset of 
1.60 kms~^ between his primary and secondary veloc- 
ities, as indicated by preliminary sine-curve fits (see Ta- 
ble 6 bv lPoDpe^l2000l where the offset he finds is similar). 
Effectively the two stars yield different center-of-mass ve- 
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locities. We found no such offset in our own measure- 
ments, but experience indicates it can sometimes appear 
wfien there is a significant mismatch between the adopted 
templates and the real stars, and if not corrected it can 
bias the semi-amplitudes when enforcing a common cen- 
ter of mass in the fit. We have thus added —1.60 kms~^ 
to Popper's secondary velocities. Finally, a third adjust- 
ment is to bring Popper's overall velocity zero point into 
agreement with ours. From trial orbital fits we found 
this required a shift of -1-0.37 kms""'^ to his velocities.^ 
Popper's corrected velocities are listed in Table [H Sep- 
arate fits to his data and ours give similar values for 
the semi-amplitudes, and yield masses that differ by less 
than twice their combined uncertainties. We therefore 
proceed to merge the two data sets below. 

A third set of ra dial velocities for CV Boo was reported 
bv lNelsonI (|2004bD . but they were obtained at lower res- 
olution, they are few in number (12), and show a much 
larger scatter than the two other data sets ('^15 kms^^) 
so they are of little use for our purposes. 

3. MODELING OF THE PHOTOMETRIC OBSERVATIONS 

The overall shape of CV Boo's light curves (Figure [2]) 
shows rather moderate proximity effects despite the sys- 
tem's relatively short period of slightly more than 20 
hours, with the curvature between the minima being 
mostly due to the deformation of the components and, 
to a smaller degree, to the mutual illumination. A num- 
ber of small-scale features are obvious to the eye that 
are possibly due to spots, other intrinsic variability, or 
even i nstrume nt al effe cts, especially in the smaller data 
sets of iNelsonI (j2004bf ). Other features described below 
are revealed through a more detailed examination, and 
introduce some complications into the analysis. 

3.1. Initial solutions without spots 

To begin with, we chose to model all the obser- 
vations together in order to obtain a baseline solu- 
tion against which to compare more complex solutions 
that attempt to account for the features mentioned 
above. We used a ve rsion of the Wilson-Devinn e y (WD) 
mode l ing program (jWilson fc DevinnevI 119711 : IWilsonI 
Il979l Il993i) with extensive modi fications as described 
m Andersen fc Claret ()2007l ) and references given 
therein. The modifications pertinent to CV Boo in- 
clude the capability to use model atmospheres (now also 
available in the distributed versions of WD), consistency 
checks between various parameter s, and the ability to use 
the downhill simplex algorithm (|Nelder fc MeadI Il965f) 
instead of differential corrections. We combined our 
own T^-b and ligh t curve with the sparser V and / light 
curves of INelsonI (|2004bD in order to improve the con- 
straint on the effective temperature ratio, and with our 
radial velocities from ii l2.3l as well as those of iPoppeil 
(j2003). Thus we solved simultaneously 3 light curves and 
4 radial-velocity curves. The parameters adjusted were 
the orbital inclination i, the secondary Tcff (flux- weighted 
mean surface temper ature) , the b andpass-specific pri- 
mary luminosity (see IWilsonlll993f ). both stellar surface 

^ While in principle the latter two adjustments (offsets) could 
be accounted for in our combined photometric and radial velocity 
solution described below by simply adding free parameters to the 
fit, limitations in the current version of the WD code do not allow 
this, so we have applied the offsets externally. 



gravitational pseudo-potentials (related to the stellar 
radii), the center-of-mass radial velocity 7, the mass ra- 
tio q = Mb /Ma, and an arbitrary phase shift. The sim- 
ilar depths of the two minima in both V and / (Fig- 
ure [2]) imply that the components must have rather sim- 
ilar temperatures, consistent with indications from spec- 
troscopy. The primary temperature was held fixed at a 
value determined from our results based on photometry 
and spectroscopy, as follows. A photometric estimate 
of the primary temperature was derived from the mean 
system temperature (i) l2.ip using approximate values for 
the radius ratio and temperature ratio from preliminary 
light-curve solutions. The result, 5755 ± 60 K, was then 
combined with the spectroscopic value of the primary 
temperature fi) l2.3|) . giving a weighted average for Star A 
of 5760 ± 150 K, which we adopt. 

The reflection albedos for both components were held 
fixed at the value 0.5, appropriate for stars with convec- 
tive envelopes, and the gravity-brightening exponents /3 
were computed internally in WD using the local value 
of Toff for each point on the stellar surf ace taking into 
accou nt m utual illuminat i on, fo llowing lAlencar fc Vaj 
(|1997t l and lAlencar et aD (|1999l ). The flux from each 
of the components is represented by NextCen atmo- 
sphere models based on the PHOENIX stellar atmo- 
sphere code (lAllard fc Hauschildtl 119951 : lAllard et all 
119971 : iHauschildt et al.lll997al |bh. The luminosity of the 
secondary is calculated by the program from its size and 
Toff- The limb-darkening coefficients for both compo- 
nents, a;A and xb, were taken from the tables bv lClaretl 
(2000), and interpolated using a bi- linear scheme for the 
current values of log 5 and Toff at each iteration. We have 
only considered the linear law here in view of the distor- 
tions in the light curve, which will tend to overwhelm 
the rather subtle effect of limb-darkening. Although we 
have no evidence of another star in the system, the pos- 
sibility of third light (^3) was explored carefully for its 
potential influence on the geometric parameters, partic- 
ularly in the solutions described below in i) l3.2i which 
include spots. Achieving convergence when solving for 
third light was found to be very difficult due to the in- 
trinsic variability and the large number of free parame- 
ters, even when considerin g multiple parameter subsets 
(jWilson fc Biermannlll976l ). We found that the solution 
was not improved, and £3 was not considered further. 
We estimate a conservative upper limit to £3 of ~1%, 
which does not produce significant changes in the geo- 
metric parameters. A circular orbit has been assumed 
in the following, based on our investigation of the eclipse 
timings in § 12.21 and the lack of any indication in the light 
curves of a displacement of the secondary minimum from 
phase 0.5. This is consistent with expectations from tidal 
theory for an orbit with such a short period (see §[6]). We 
have also assumed here tidal forces have synchronized the 
components' rotation with the orbital motion of CV Boo. 
In these initial calculations we applied both least-squares 
differential corrections and/or the simplex method be- 
tween successive iterations. The limb-darkening coeffi- 
cients, normalization magnitudes, surface gravities, and 
individual velocity amplitudes were all updated between 
consecutive runs to correspond to the solution from the 
previous iteration. 

The solutions are shown in Figures [2l [3] and IH with the 
corresponding residuals. The residuals from the radial 
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Fig. 3. — Primary (x) and secondary (+) radial velocity mea- 
surements collected at CfA (Table [3j along with the theoretical 
curves obtained with WD and no spots (top panel). Velocities are 
shifted so that the center-of-mass velocity 7 is at zero (dashed line). 
The large deviation from Keplerian motion in the predicted veloc- 
ity around both conj unction s i s due to the Rossiter-McLaughli n 
effect (see ISchlesingerlll909ai rbl: IRossite^lT92l; iMcLauehlinI IT92I) . 
caused by partial eclipses of the rotating stellar surfaces (an ef- 
fect built into the WD model). The O — C residuals are shown at 
the bottom. The standard deviation of the unweighted residuals 
is (Trv = 6.09 kms~^ for both components and is shown in the 
upper left corner. The reduced values were 1.000 and 0.998, 
respectively. 

velocity fits match the quality of the observations. The 
photometry is reasonably well represented on average by 
the theoretical curves, but the residuals for our F-band 
observations show an rms scatter of 0.0196 mag that is 
much larger than the mean internal error (~0.013 mag). 
This immediately suggests there may be unmodeled ef- 
fects. The intrinsic errors of Nelson's observations were 
not reported in the original publication. 

3.1.1. Study of the light-curve residuals 

Part of the extra scatter is no doubt due to features in 
the light curve alluded to earlier that are seen in Figure[21 
such as changes in the light level from night to night at 
the same orbital phase (e.g., near phase 0.1), or other 
short-term deviations (e.g., near phase 0.8), both in our 
data and in Nelson's. We investigated the residuals of 
our more numerous y-band photometry further to search 
for periodic signals that might additionally contribute to 
the excess scatter. Our initial explorati on of possible sig- 
nals n ear the orbital period using the iLafler fc KinmanI 
(|1965f ) method revealed several similar periodicities that 
appear significant. We then extended the search to a 
much wider range of frequencies by computing the Lomb- 
Scargle periodogram, and found other signals. This is 
shown in the top panel of Figure [H] We refer to this 
as the "dirty" power spectrum, since it is affected by 
the particular time sampling of the observations (win- 
dow function). The highest peak corresponds to a pe- 
riod of ~0.837 days, which is shorter than the orbital 
period of 0.846993420 days. The two next highest peaks 
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Fi g. 4. — Prima ry (x) and secondary (-|-) radial velocity data 
from lPoppeil 1 I2OOOI) . with the same limits on the vertical axes as in 
Figure |3] The standard deviations of the unweighted residuals are 
""rv.A = 3.14 kms~^ and o"rv,B = 2.14 kms~^ for the primary and 
secondary velocities, respectively (shown in the upper left corner 
of the lower panels). The reduced values were 1.005 and 1.002. 

(indicated with arrows) turn out to be 1-day aliases. To 
illustrate this, we have applied t h e CL EAN algorithm 
as implemented by iRoberts et al.l (|1987D to remove the 
effects of the window function. The second panel of Fig- 
ure [5] shows that only the 0.837-day peak survives this 
process, suggesting it is a real signal. In the third panel 
an enlargement of the dirty power spectrum in the vicin- 
ity of the mai n peak reveals fine struct ure that was also 
seen with the ILafler fc KinmanI ()1965f ) method. How- 
ever, none of these peaks agree with the frequency corre- 
sponding to the orbital period, which is represented for 
reference with a dotted line. The two main sidelobes in- 
dicated with arrows are 1-year aliases of the main peak. 
Once again they disappear after application of CLEAN, 
as seen in the bottom panel, supporting the reality of 
the remaining signal. The statistical significance of this 
signal was estimated by numerical simulation. We gener- 
ated 100,000 artificial data sets using the actual times of 
observation and the variance of the original residuals as- 
suming a Gaussian distribution of errors, and computed 
the Lomb-Scargle power spectrum for these data sets over 
the same frequency interval considered above. We then 
selecting the highest peak in each case. None of them 
came close to the height of the peak we see in the real 
data, indicating a false alarm probability smaller than 
10-5. 

The precise frequency of this signal was measured in 
the CLEANed spectrum, and its uncertainty was es- 
timated from the half width at half maximum of the 
peak. The corresponding period is 0.83748 ± 0.00052 
days, which is different from the orbital period at the 
18(7 level. A plot of the photometric residuals folded 
with this period is shown in Figure [H and indicates a 
peak-to-peak amplitude of about 0.04-0.05 mag. 

The Kimpel Observatory data cover two observing sea- 
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Fig. 5. — (a) Lomb-Scargle power spectrum of the residuals 
of our V-band Kimpel Observatory observations of CV Boo from 
the no-spot solution described in the text. The arrows indicate 
1-day aliases of the central peak; (b) CLEANed power spectrum 
of the same measurements using the algorithm of .Roberts et al. 
II1987I) to remove the effects of the time sampling; (c) Enlargement 
of panel (a), with the 1-year aliases of the main peak indicated 
with arrows. The dotted line represents the orbital frequency; (d) 
Enlargement of panel (b). The period corresponding to the signal 
is 0.83748 ± 0.00052 days. 
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Fig. 6. — Residuals of our V-band Kimpel Observatory observa- 
tions of CV Boo shown as a function of phase, using the period of 
0.83748 days inferred from the power spectrum analysis (see Fig- 
ure|5}. The time origin has been set arbitrarily to HJD 2,450,000. 



sons. Separate Lomb-Scargle power spectra show that 
the same signal is present in both seasons, along with 
the 1-day aliases, indicating the phenomenon is persis- 
tent from one year to the next. It is not seen as clearly, 
however, in t he residu als from our baseline fit of the ob- 
servations of iNelsonI ([2004b ) , which are much sparser 
than ours (and span only 71 days instead of 549 days). 
His T^-band data show a hint of the main peak and its 1- 
day aliases, but not the /-band data, which have a larger 
scatter. 

We carried out a similar power spectrum analysis of 
the (comp— ck) differential magnitudes from Kimpel Ob- 
servatory, to explore the possibility that either the com- 
parison or the check star might be the source of this 
variation. No significant periodicity was seen. Thus, 
the phenomenon is intrinsic to CV Boo. Possible expla- 
nations for this variation include stellar pulsation, and 
star spots on one or both components. Neither of the 
stars appear to be in an evolutionary state that favors 
pulsational instability. For example, the absolute dimen- 
sions derived below place both components well outside 
the Cepheid or 6 Set instability strips in the H-R dia- 
gram indicated by |K]^rgaard£tin] (Illll). On the other 
hand, CV Boo is a known str ong X-ray source detected 
by ROSAT (jVoges et al.ll2000D . with an X-ray luminosity 
of logix = 30.65^ that is some 4000 times stronger than 
the Sun. In terms of its bolometric luminosity CV Boo 
has logLx/iboi = —3.39, which is near the high end 
for active binaries. Thus it seems likely that the under- 
lying reason for the 0.837-day periodicity is related to 
spots on the surface of one or both components, and we 
proceed under this assumption. It is interesting to note 
that these features on CV Boo seem to have lasted for 
an unusually long time (1.5 years in our case), at least 
compared to sunspots, although even more extreme ex- 
amples have been documented in the literature. One i s 
the well-known active binary HR 1099 (|Vogt et al.lll999( l. 
with surface features persisting for at least 11 years. 

An important implication of this spot hypothesis is 
that the component having spots would appear to be 
rotating slightly more rapidly than synchronously with 
the motion in the circular orbit, which is unexpected for 
such a short-period binary. We discuss this in more detail 
below. 

3.2. Solutions with spots 

An accurate measurement of the v sin i values for both 
components would allow for a direct test of our hypoth- 
esis of non-synchronous rotation, and could even distin- 
guish which of the stars is the culprit (or if both are). 
Unfortunately, however, the quality of our spectroscopic 
material is insufficient for that purpose. In principle, 
modern light-curve models such as WD enable the user 
to solve for various parameters that describe the spots. 
However, with only photometric data at our disposal for 
CV Boo, and most of it in a single passband, it is es- 
sentially impossible to tell which star has the spots, or 
whether both components have them. This is a well- 
known difficulty in light-curve modeling. Other inver- 
sion techniques such as Doppler imaging are much bet- 

Z/x is in units of erg , and was determined from the ROSAT 
count rates and hardness ratio s, the distance estima te in §|4] and 
the energy conversion factor of lFleming et al.l HlMB ). 
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ter suited to mapping surface inhomogeneities, although 
even they are not without their Hmitations. Moreover, 
even if we knew which star has the spots, the determina- 
tion of their parameters from Hght curves alone is a no- 
toriously ill-posed problem, on which there is abundant 
literature discussing issues of indeterminacy and non- 
uniq ueness in the pre sence of limited data quality (see, 
e.g.. [Ekei1[l996Lfl999l and numerous references therein). 
Having photometry in multiple passbands may aleviate 
the problem somewhat, but it doesn't solve it and strong 
degeneracies are likely to remain with other subtle effects 
in the light-curves. Therefore, while we cannot hope to 
obtain an accurate picture of the distribution of any sur- 
face features here, the consequences of spots on the light 
curve are fairly clear in CV Boo (to the extent that our 
hypothesis is true), and we make an effort in the fol- 
lowing to at least remove some of those distortions and 
study their influence on the geometric parameters of the 
system, which are of more immediate interest. 

In order to permit the numerical treatment of surface 
features in this case, we introduced modifications in the 
WD code to allow for a precise tracking of the spot po- 
sition at a period different than the orbital one. In this 
scheme, we specify the spot properties at a certain Julian 
date and, through the specified intrinsic rotation rate, 
the code keeps track of the spot motion, with its lon- 
gitude following the component's rotation, and its co- 
latitude, size, and effective temperature remaining oth- 
erwise constant. In view of the ambiguities mentioned 
above regarding the location of spots in CV Boo, and 
our inability to tell if there might even be multiple spots 
on one or both stars, we have taken our light-curve fit 
from § 13.11 as our starting point and investigated the fol- 
lowing three simple cases separately: (a) a single spot 
on the primary; (b) a single spot on the secondary; and 
(c) one spot on each component. More complex con- 
figurations become increasingly difficult to study due to 
convergence problems in the solutions, and it is not clear 
they are justified with the data available. 

The influence of spots on the radial velocity curves is 
very small compared to our errors, so that those data 
are not of very helpful for studying surface features. We 
use only our more extensive V-band photometry in the 
study of these three cases, although the solutions were 
checked using Nelson's V and / light curves. As the pho- 
tometric coverage does not necessarily overlap with the 
radial velocity coverage, we have adopted for the spotted 
cases the mass ratio q obtained from a no-spot fit similar 
to that in that assumes asynchronous rotation (see 
below) , and held it fixed. For lack of other physical con- 
straints, and given that the stars are quite similar in all 
their properties, we assumed that both components ro- 
tate slightly super-synchronously, at a rate given by the 
ratio between the orbital period and the residual period 
found in the previous section, equal to 1.0114. 

Cases (a) and (b) converged quite rapidly to similar 
configurations, in which the co-latitude, size, and the ef- 
fective temperatures of the spots are comparable, while 
their longitudes are such that the spots present always 
the same position relative to the center of the spotted 
component and the observer. This can be seen in TablelHl 
where the longitudes of the spots in cases (a) and (b) are 
separated by nearly 180°. Another similarity is that both 
spots cover the components' polar regions. Although at- 
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Fig. 7. — Difference between the light curves with spots and 
those without spots (shown in Figure 01, for the first cycle of our 
observations. The thin horizontal line shows the normalization 
level of the theoretical light curves (phase 0.25, first orbiting cycle). 



tempted, no solution could be obtained for "hot" spots 
(i.e., with temperature factors larger than unity; see be- 
low). 

Solution (c) with one spot on each component did not 
converge as easily. When the parameters of both spots 
were left free to be adjusted, one of the spots (usually 
the one on the primary) tended to become very small 
and cold, with the temperature factor Tfactor (ratio be- 
tween the spot temperature and the photospheric tem- 
perature) becoming smaller than allowed by the NextGen 
atmosphere tables we used. The solution we present was 
achieved by first adjusting some of the spot parameters 
while holding others fixed, and then alternating and it- 
erating until convergence. 

The maximum amplitude of the influence of the spots 
on the light curve is '--^0. 08 mag, and occurs for the one- 
spot solutions, as shown in Figure [Tj This figure corre- 
sponds to the first orbital cycle of our observations and, 
since the spots follow the components' non-synchronous 
rotation, the dips change place at each orbiting cy- 
cle. The two-spot solution of case (c) gives a slightly 
smaller peak-to-peak amplitude (^0.06 mag) that seems 
marginally larger than indicated in Figure [HI suggesting 
that perhaps a more complex spot configuration may be 
needed. 

We report in Table [6] the model parameters we ob- 
tain for the solution with no spots and for cases (a), 
(b), and (c), together with the radii of the components 
in terms of the orbital separation. For the reasons de- 
scribed above, the solution without spots was performed 
by solving simultaneously three light curves and 4 radial 
velocity curves, whereas the spotted fits are based only 
on our T^-band light curve. The main difference in the 
parameter values is seen in the inclination angle, which is 
approximately one degree higher for the solution without 
spots. Other parameters such as the secondary effective 
temperature and the sizes of the components tend to dif- 
fer less between the spotted and unspotted solutions. 

Figure [8] gives a representation of the spot configura- 
tion resulting from case (c) with the components' size 
and separation rendered to scale, and seen from the ob- 
server's viewpoint at six different orbital phases. The 
stars are well detached fror n the corre sponding Roche 
lobes, with fill-out factors (jMochnackil [19841 that are 
0.7693 and 0.7501 for the primary and secondary, respec- 
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Fig. 8. — Representation of the components of CV Boo at differ- 
ent orbital phases as indicated on the left, shown to scale with their 
high-latitude spots as modeled here. These spots resulting from 
our solution (c) are positioned in such a way that they practically 
avoid being eclipsed. This is an indication of a rather sinusoidal 
behavior of the disturbing phenomenon causing periodic variations 
in the residuals of the unspotted solution (see text). 

tively. We noted above that our fits yield polar spots, as 
has often been found (also from Doppler imaging tech- 
niques) for other active binaries such as the RS CVn sys- 
tems. There is considerable theoretical support for this 
preference for high-latitude surface features in rapidly- 
rotati ng active systems ( s ee, e.g.. iSchiissler fc Solankil 
Il992t iGranzer et al.l l200a ll^ik et all l2007f) . A curious 
result from our fits is that the spots happen to be po- 
sitioned so as to avoid eclipses, although the reality of 
this configuration is difficult to assess. It is nevertheless 
an indication that the phenomenon responsible for the 
periodic behavior of the residuals in the unspotted solu- 
tion does not lead to strong discontinuities, such as those 
resulting from the eclipses. 

Although the rms residual is marginally smaller for 
the solution obtained in case (c), as indicated in Table [SI 
this fit as well as the other two spotted solutions are vi- 
sually indistinguishable from the solution without spots 
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Fig. 9. — O — C residuals from the solution with one spot on 
each component, based on the fit to our V-band light curve. The 
model used to compute the residuals for the V and / light curves 
of Nelson (2004b) in the top panels is based on the same light- 
curve parameters as our V-band fit, except for the magnitude at 
quadrature and the wavelength of the observations. The standard 
deviation of a single observation for each residual curve is shown 
in the upper left corner of each panel. 




3 - : I — 

O.SO 0.B2 0.B4 0.S6 CBS 0.90 

Trial period (days) 

Fig. 10. — Th e variance versus trial period in days (following 
ILafler fc Kinmaiil lT965) for the Kimpel Observatory V-band resid- 
uals of the no-spot solution (top curve, corresponding to the fit 
shown in Figure [2J, and for the solution with two spots (bottom 
curve; see text). The dotted line marks the most significant period 
found for the O—C of the solution with no spots, while the dashed 
line indicates the orbital period. 

shown in Figure [2l The residuals for case (c) are dis- 
played in Figure [3] for our y-ban d light curve as well 
as for the V and I hght curves of lNelsonI ()2004bD . The 
patterns clearly visible in these O — C diagrams are not 
very different from those in Figure [21 which may give 
the impression that not much progress has been made.^ 
They certainly indicate that there are still features of the 
brightness variation that are not completely or correctly 
modeled, possibly due to a more complex spot configu- 
ration than we have assumed, or even some combination 
of spots and multi- modal pulsations. Problems of an in- 
strumental nature in the photometry cannot entirely be 
ruled out either. However, what is not immediately obvi- 
ous to the eye is that no significant periodicities that we 
can detect remain in these residuals. This is illustrated 
in Figure I10[ in which the top curve shows our Lafler- 
Kinman period study of the Kimpel Observatory 1/^-band 
residuals from the no-spot solution, and the lower curve 
shows the same study for the residuals from case (c). 
Note the common vertical axis for both sets of residuals, 
indicating the improvement in the overall variance. 

4. ABSOLUTE DIMENSIONS AND PHYSICAL PROPERTIES 



Note, however, that those patterns are most obvious in Nelson's 
data, which do not actually enter into the final solution adopted 
in §|4] 
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Examination of Table [6] shows that key geometric pa- 
rameters such as the relative radii (rA,voi, ''b,voi) vary by 
as much as 3-4% between the three spotted solutions, 
with solution (c) generally giving intermediate results. 
As indicated earlier, this is the fit that provides formally 
the smallest rms residual, although the difference com- 
pared to the other two spotted solutions is marginal. In 
all four solutions the mean light ratio outside of eclipse 
accounting for spots, (^b/^a)^: is quite similar to the 
spectroscopically determined value of 0.73 ± 0.04 (i) l2.3p . 
From the effective temperatures of CV Boo A and B the 
convective turnover tim e for both stars is estimated to 
be ~25 days, following (|1994[ ). The Rossby num- 
ber (ratio between the rotation period and the convec- 
tive turnover time) is then Rq « 0.033, which places 
both components in the regime where stars usually dis- 
play significant light variations due to spots (see iHall 
11994 Figure 6). On the basis of the above we adopt fit 
(c) with one spot on each component as the best com- 
promise for CV Boo, but we reiterate that this model 
is still probably only a crude approximation to the true 
spot configuration in the system, assuming that spots 
are the underlying reason for the periodic signal found 
in the light-curve residuals. For calculating the absolute 
dimensions of the two stars we have chosen to use more 
conservative uncertainties than the formal errors listed in 
Table [6l to account for the spread among the three spot- 
ted solutions given the uncertainties in the modeling: we 
have combined the internal errors quadratically with half 
of the maximum range in each parameter. The values 
adopted are i = 86?24 ± 0?33, a = 4.748 ± 0.019 R©, 
q = 0.9378 ± 0.0070, ta.voI = 0.2658 ± 0.0047, and 
'"B,voi = 0.2470 ± 0.0048, and are based only on the Kim- 
pel Observatory measurements. The final results are pre- 
sented in Table [71 where the uncertainties were obtained 
by propagating all observational errors in the usual way. 

The stars in CV Boo depart somewhat from the spher- 
ical shape due to tidal and rotational distortions. The 
relative difference between the polar radius and the ra- 
dius toward the inner Lagrangian point is 5.5% for the 
primary and 4.8% for the secondary. The system is nev- 
ertheless well detached: the sizes of the stars represent 
fractions of 70% and 66% of their respective mean Roche 
lobe sizes. The temperature for the secondary from the 
light-curve solution is in excellent agreement with the 
spectroscopic value fii l2.3p . 

Included in Table [7] are the predicted projected rota- 
tional velocities (tiasync sini) computed with the adopted 
rotation period for the stars (Prot = 0.83748 days = 
Porb/1-0114; see § I3.1.1|) . as well as the synchronous 
values (iigync sinz), for reference. These may be com- 
pared with the measured v sin i values from spectroscopy 
(ij |2.3[) . The stellar radii used for these calculations are 
those presented to the observer at quadrature (which are 
2.7% and 2.4% larger than the volume radii; see Table[6|), 
since that is the phase at which the spectroscopic obser- 
vations are concentrated. As a proxy for the radius at 
quadrature we use the average of rpoint and rback- 

Finally, for computing the absolute visual magnitude 
My and distance we have relied on the apparent V mag- 
nitude listed in the Tycho Catalog, and ignored extinc- 
tion. CV Boo was not o bserved by the Hipparcos mission 
(jPerrvman et al.|[l997[ ). so no direct parallax measure- 
ment is available. 
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Fig. 11. — Absolute dimensions for CV Boo com pared with evo- 
lutionary models from the series bv IYi et al.lll200ll '). The error bars 
for log g are smaller than the size of the symbols. Mass tracks for 
the exact masses we measure are indicated with solid curves for 
the best-fitting metallicity of Z = 0.01955 (where Zq = 0.01812 
for these models). Solar metallicity tracks are shown for reference 
(dotted curves). The isochrone producing the best simultaneous 
match to both components is shown with the short-dash line, and 
corresponds to an age of 9.0 Gyr. The long-dash lines represent 
small sections of the two isochrones corresponding to the maximum 
and minimum age allowed by the errors (9.0 ± 1.8 Gyr). The un- 
certainty in the location of the mass tracks is indicated with the 
error bar below the tracks for the primary, and is much smaller 
than the temperature uncertainty. 

5. COMPARISON WITH STELLAR EVOLUTION THEORY 

In this section we compare the absolute dimensions 
of CV Boo with current stellar evolutio n models from 
the Yonsei-Yale series bv lYi et al.l ()2001[ ). incorporating 
an updated pres cription for convective c ore overshooting 
as described bv iDemargue et al] (120041 ). These models 
adopt a mixing length parameter of auL = 1.7432, cal- 
ibrated against the Sun. In Figure [TT] we show evolu- 
tionary tracks computed for the exact masses we derive 
for each star (see Table [7]) , for a heavy-element abun- 
dance equal to that of the Sun (which is Zq — 0.01812 
in these models; dotted lines). The uncertainty in the 
location of the tracks that comes from our mass errors 
is indicated with the error bar in the lower left. The 
tracks show excellent agreement with the observations, 
suggesting the composition is near solar. The measured 
temperature difference between the components is quite 
close to what the models predict. A marginally better 
match is achieved with a slightly higher abundance of 
Z = 0.01955 (corresponding to [Fe/H] = +0.04, assum- 
ing no enhancement of the a elements), shown as solid 
lines in the figure. The models indicate the primary is 
beginning its shell hydrogen-burning phase, and the sec- 
ondary is near the end of its main-sequence phase. The 
age that best fits both components in this log g-T^g di- 
agram is 9.0 ±1.8 Gyr, and the corresponding isochrone 
is shown as a dashed line. 

We have al so considered a second set of models, from 
the series by IClaretl ()2004|) . The physics in these cal- 
culations is similar though not exactly the same as the 
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Fig. 12. — Absolute dimensions fo r CV Boo co mpared with 
evolutionary models from the series bv lClaretl I I2004I) with a value 
of the mixing-length parameter of qml = 1.68. Mass tracks for the 
exact masses we measure and for solar composition (Z = 0.020, 
X = 0.70, in these models) are indicated with dotted curves. The 
solid curves giving a somewhat better fit correspond to models with 
a slightly higher metallicity oi Z = 0.023. An isochrone for log age 
= 9.99 is shown for reference (dashed line). The uncertainty in the 
location of the mass tracks is indicated with the error bar below the 
tracks for the primary, and is much smaller than the temperature 
uncertainty. 

previous ones. For example, the solar composition in this 
case is taken to be Zq — 0.020, and the mixing length 
parameter that best reproduces the observed properties 
of the Sun is q;ml = 1.68. The comparison with the 
observations for CV Boo is shown in Figure [T^l Al- 
though the Claret models match the measured proper- 
ties very well, we find as with the Yonsei-Yale models 
that a shghtly higher metaUicity (Z = 0.0225, or [Fe/H] 
= -1-0.05) provides an even better fit. This is shown by 
the solid lines in Figure [T^l The age of the system from 
these calculations is 9.8 Gyr, consistent with the previ- 
ous estimate. Experiments changing the mixing length 
parameter show the sensitivity of the best-fit composi- 
tion to ttML • In Figure [T3] we compare the observations 
with tracks computed for a lower value of oml = 1-50, 
which has the effect of yielding lower temperature pre- 
dictions. Solar-metallicity models are indicated with the 
dotted lines. In this case we find that the best-fit metal- 
hcity {Z = 0.0185, or [Fe/H] = -0.03) is slightly lower 
than solar (solid lines). 

The preceding comparisons may give the impression 
that the observations for CV Boo are very well matched 
by the predictions from theory, that stellar physics is well 
understood, and that therefore there is no reason for con- 
cern. However, a more careful examination indicates that 
this is not necessarily true. Of the three basic parameters 
typically determined in eclipsing binaries (M, i?, Tcs), 
the temperature is usually the weakest since it often relies 
on external calibrations. Figure [TJ] displays the measure- 
ments for CV Boo in a different diagram, the mass-radius 
plane, along with isochrones from the Yonsei-Yale series 
for the same two metallicities discussed in Figure [TTJ No 
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Fig. 13. — Same as Figure [TJ] but for a mixing-length parameter 
of aML = 1.50. Mass tracks for the measured masses and for 
solar composition {Z = 0.020, X = 0.70) are indicated with dotted 
curves. The somewhat better-fitting solid curves correspond to 
a slightly lower metallicity of Z = 0.0185 in this case, showing 
the influence of the qml parameter in the determination of the 
composition of CV Boo. 

single model matches both components within the errors, 
and the secondary appears nominally older than the pri- 
mary, the difference in age bei ng ^25%. This is the same 
phenomenon pointed out by iPoppeH ()1997[ ) for several 
other systems including FL Lyr, RT And, UV Psc, and 
a Cen. Another way of interpreting this is that the secon- 
daries in all these binaries are too large for their masses, 
compared to theory or compared to the primaries. For 
CV Boo the offset in the secondary radius is ~10%, 
which represents a very significant 5(t deviation. Sim- 
ilar radius discrepancies have been described recently by 
others (e.g. , jCla us cn et al. 19 99a: Torres & Ribas 200j; 
iRibad [20031: lL6pez-Morales fc RibasI 120051: ITorresll20oC 
although ear ly in dications go as far back as the work of 
iHoxid (|l973[ ) and iLacvl ()1977[ ). The prevailing explana- 
tion seems to be that the enlarged radii of the secon- 
daries, which are typically well under a solar mass, are 
caused by strong magnetic fields and/or spots commonly 
associated with chromospheric ac t ivity i n these systems 
(see e.g., iMullan fc MacDonaldl 120011: iChabrier et~aLl 
120071 for the theoretical context). The signs of activity in 
CV Boo are fairly obvious (spottedness. X-ray emission), 
and are no doubt associated with the rapid rotation of 
the components. 

6. COMPARISON WITH TIDAL THEORY 

The predictions of tidal theory were compared with 
the observations by computing the time of circular- 
ization and synchronization fo r C V Boo us ing t he ra- 
diativ e damping formalism of IZahnI ()1977| ) and IZahnI 
(1 19891). as well as the hydrodynamical mechanism of 
iTassoul fc Tassoul ()1997D . and references t herein. The 
proced ure fo llows closely that describ ed bv lClaret et al.l 
(jl995[ ) and IClaret fc Cullhal (|1997[ ). Both theories 
predict that synchronization and circularization are 
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Fig. 14. — Mass-radius diagram for CV Boo, showing the mea- 
surements against isochrones from the Yo nse i-Yale series for the 
same two metallicitics displayed in Figure [TTI Ages are indicated 
along the top. 

achieved very quickly in this system by virtue of the short 
orbital period, at an age of merely 157 Myr (logt = 8.197, 
or less than 2% of the evolutionary age). The fact that we 
measure the orbit to be circular is therefore not surpris- 
ing. On the other hand, the evidence from our photomet- 
ric observations (§ l3.1.ip suggesting the rotation may be 
slightly super-synchronous for at least one of the compo- 
nents is more interesting, as it is not predicted by theory. 
Given the nature of the system, an activity-related ex- 
planation to this discrepancy is certainly possible. More 
precise measurements of the projected rotational veloci- 
ties V sin i for the components would be very helpful. 

7. DISCUSSION AND CONCLUSIONS 

Despite the system's intrinsic variability, the absolute 
dimensions for the components of CV Boo have now been 
established quite precisely. The relative errors are bet- 
ter than 1.3% in the masses and 2% in the radii. The 
object can now be counted among the group of eclipsing 
binaries with well-known parameters. Under different 
circumstances the large number and high quality of the 
photometric observations we have collected might have 
permitted a more detailed study of the limb darkening 
laws and a comparison with theoretically predicted co- 
efficients, but this possibility was thwarted here by the 
intrinsic variability. This phenomenon is not itself with- 
out interest. If interpreted as due to the presence of 
spots, as we have done here, it implies that at least one 
of the stars is rotating about 1% more rapidly than the 
synchronous rate, a result that was unexpected for a 
close but well detached system such as this. We con- 
clude that our current understanding of tidal evolution 
is still incomplete, or that other processes are at play 
in this system that theory does not account for. One 
interesting possibility is differential rotation. The inter- 
pretation of measurements of the rotation period of a 
star made by photometric means, as we have implicitly 
done here, usually relies on the assumption of solid-body 
rotation. More often than not, spots are located at inter- 
mediate latitudes rather than on the equator, or at high 



latitudes in more active stars, and differential rotation 
is such that the stellar surface revolves more slowly at 
higher latitudes, at least in the Sun. This will tend to 
bias photometric rotation measurements towards longer 
periods, if differential rotation is significant enough. In 
CV Boo we see the opposite: the period is shorter than 
the equatorial rate, assuming that synchronization holds. 
Thus differential rotation can only explain the signal we 
have detected if it is "anti-solar" , with the polar regions 
rotating more rapidly. A handful of stars do indeed show 
evidence of weak anti-solar differen tial rotation (e.g. , 
IL Hya, HP 31933, a Gem. UZ Lib: IWeber et al.l 120031: 
Strassmeier et all 120031 : iKSvari et all l2007t IVida et all 
2007( ). Thev all happen to be very active (some of them 
with high-latitude spots, as in CV Boo), although they 
tend to be giants or subgiants rather than dwarfs. It 
is thought that this phe nomenon may result from fas t 
meridional flows (see, e.g.. lKitchatinov fc Riidigei1l2004f) . 
Further progress in understanding the rotation of the 
CV Boo components could be made with additional dif- 
ferential photometric observations in several passbands, 
along with simultaneous high-resolution, high signal-to- 
noise ratio spectroscopy over a full orbital cycle. 

Another significant discrepancy we find with theory 
is in the radius of the secondary, which appears to be 
some 10% too large compared with predictions from stel- 
lar evolution models. This difference is in the same di- 
rection as seen for a number of other low-mass eclips- 
ing binaries, as mentioned in §[5l In those cases one of 
the explanations most often proposed is that the strong 
magnetic fields associated with activity (which is com- 
mon in rapidly-rotating K and M dwarfs in close bina- 
ries) tend to inhibit convective motions, and the struc- 
ture of the star adjusts by increasing its size to allow the 
surface to radiate the same amount of energy. At the 
same time, the effective temperature tends to decrease 
in order to preserve the total luminosity. Spot cover- 
age can produce similar effects. Theoretical and obser- 
vational evidence for the conservation of the luminosity 
in these systems has been presented by Delfosse et ai] 
' 200Q) . iMuUan & MacDonald (2001), Torres & Ribai 
2002 ). lRibasl (|2006l) . ITorres et al. (2006), Cha brier et all 
2007t i. and others (see also [ Morales et al.. 20081 ). We do 



not see any obvious discrepancy in the temperature of 
CV Boo B compared to models, although our uncertain- 
ties are large enough that the effect may be masked. 

If we restrict ourselves to well studied double-lined 
eclipsing binaries in which the mass and radius deter- 
minations are the most reliable, deviations from theory 
such as those described above have usually been seen in 
stars that are considerably less massive than the Sun, 
which have de ep convective envelopes. However, the re- 
cent study bv [Torres et al.l ()2006[ ) pointed out that the 
problem is not confined to the lower mass stars, but 
extends to active objects approaching 1 M0, such as 
V1061 Cyg Ab, with M = 0.93 Mq. CV Boo B has 
an even larger mass of 0.968 Mq, and also appears to be 
oversized. Similarly with the virtually identical active 
star FL Lyr B {M = 0.960 Mq). The convective en- 
velopes of these objects are considerably thinner than in 
K and M dwarfs and represent only a few percent of the 
total mass, yet they appear sufficient for magnetic fields 
to take hold and alter the global properties of the star, 
if that is the cause of the discrepancies. These examples 
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show once again that our understanding of steUar evolu- 
tion theory is incomplete, even for stars near the mass of 
the Sun. 
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TABLE 1 

Differential V-band 
measurements of cv boo. 



HJD-2,400,000 


Phase 


AF 


52250.99816 


0.35467 


+0.452 


52250.99907 


0.35575 


+0.416 


52251.00000 


0.35685 


+0.445 


52251.00091 


0.35792 


+0.508 


52251.00182 


0.35900 


+0.433 



Note. — Table [T] is available in its 
entirety in the electronic edition of the 
Astronomical Journal. A portion is 
shown here for guidance regarding its 
form and contents. 



TABLE 2 

Photometric indices and inferred mean effective 
temperature of cv boo. 



Photometric Index Value T^s (K) Ref. 



Johnson V 




10.75 ± 0.10 




1 


Johnson B — V 




0.73 ± 0.11 


5417 ± 350 


1 


Tycho-2 Bt - Vt 




0.82 ± 0.13 


5448 ± 329 


1 


Johnson/2MASS V- 


■J 


1.18 ± 0.10 


5693 ± 103 


1,2 


Johnson/2MASS V- 


H 


1.47 ± 0.10 


5666 ± 155 


1,2 


Johnson/2MASS V- 


Ks.... 


1.55 ± 0.10 


5692 ± 157 


1,2 


Tycho-2/2MASS Vr 


-Ks.. 


1.629 ± 0.081 


5679 ± 129 


1,2 


Sloan g — r 




0.473 ± 0.002 


5760 ± 100 


3 



Note. — References: a) IH0g et ani|2000h : (2) ICutri et al] (12003) : 

(3) Sloan Digital Sky Survey data. 
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TABLE 3 

New radial velocity measurements of CV Boo. 



HJD 
-2 440 000 


Phase 


Star A 

km / s 


(CM3) 
A 


Star B 

km/s 


(O-C) 
B 


HJD 

-2 440 000 


Phase 


Star A 

km/s 


(O-C) 

a 


Star B 

km / s 


(O-C) 

B 


4S4nS SSS1 


0.1792 


— 127 


53 


—2 


12 


+ 144 


05 


+13.68 


52805 7304 





2974 


— 133 


61 


— 1 


12 


+134 


43 


—3 


53 




0.6453 


+113 


28 


+6 


33 


— 121 


89 


—4.64 







6180 


+79 


11 


— 12 


00 


—96 


25 


+4 


12 




0.9148 


+70 


53 


+1 


73 


-80 


45 


—4.19 







7833 


+141 


29 


+8 


91 


— 158 


46 


— 14 


00 




S5'?0 


+101 


70 


—6 


29 


— 121 


62 


—3.39 







3715 


— 116 


79 


— 16 


14 


+112 


11 


+8 


13 




U. OOZrU 


—96 


29 


—2 


48 


+97 


00 


+0.31 


'^9^^'^n 7704 





8608 


+104 


81 


+0 


97 


— 115 


59 


— 1 


81 


'JinJC'-L.OiJiiJ 




+115 


23 


+4 


57 


— 126 


58 


—5.38 


^98Q4 6202 





2449 


— 148 


22 


—9 


80 


+145 


71 


+1 


37 


^23Q^ 77Sn 


0.2884 


— 126 


96 


+7 


58 


+131 


35 


—8.81 


'^'^ni 1 c\^A'\ 





7124 


+138 


93 


+7 


38 


— 147 


38 


—3 


84 




0.7395 


+133 


04 


— 1 


97 


— 156 


34 


—9.07 


5Sni 7 OfS'^fi 





8097 


+130 


12 


+4 


16 


— 144 


86 


—7 


31 




88S2 


+80 


62 


—9 


96 


— 102 


38 


—2.82 


tJOUOU. UIO^ 





2214 


— 130 


15 


+6 


19 


+149 


26 


+7 


14 




0.7022 


+119 


97 


—9 


29 


— 137 


99 


+3.09 


'\'^C\A^ {\9(\'\ 

OOyj'-ttj.Kj^KjO 





8143 


+110 


06 


— 14 


40 


— 137 


98 


—2 


03 


^24S1 72SQ 


0.7601 


+131 


60 


—3 


44 


— 147 


89 


—0.59 







3184 


— 133 


81 


-7 


71 


+138 


58 


+7 


46 


^2^S7 fim ^ 


0.7319 


+145 


12 


+10 


69 


— 153 


73 


— 7.10 


'^'^079 091 d 





6931 


+120 


32 


-6 


41 


— 138 


36 


+0 


02 




0.7327 


+136 


61 


+2 


10 


— 143 


34 


+3.38 







1890 


— 128 


08 


+0 


67 


+140 


61 


+6 


65 




0.2150 


-139 


26 


—3 


99 


+149 


89 


+8.92 







2227 


— 148 


22 


— 11 


69 


+142 


32 


— 


01 




0.1853 


— 126 


70 


+0 


86 


+125 


61 


—7 06 


"^•^l 9^^ 9 





2825 


— 135 


40 


+0 


26 


+136 


82 


—4 


54 




0.3016 


— 136 


00 


—4 


61 


+134 


21 


—2.57 


531 31 8024 





2733 


— 148 


17 


— 11 


15 


+142 


70 


— 


12 




0.7903 


+137 


79 


+6 


78 


— 152 


45 


—9.46 


00i-00.0^0\J 





6648 


+113 


70 


—2 


64 


— 125 


78 


+1 


49 


^271 2 (1278 


0.6677 


+117 


71 


+0 


12 


— 132 


03 


—3.42 


SSI S4 snos 





8134 


+123 


52 


— 1 


25 


— 134 


88 


+1 


40 


^271 8 Q482 


8383 


+103 


39 


— 11 


69 


— 126 


31 


—0.46 


Q1 o'i 





7362 


+135 


62 


+0 


82 


— 159 


54 


— 12 


50 


•^9790 QQIT^ 


0.2493 


— 142 


32 


—3 


84 


+135 


44 


—8.97 


00i-0\J.0\J i o 





7959 


+124 


95 


—4 


81 


— 145 


87 


—4 


23 


'^9791 S8fSfi 


0.3075 


-131 


97 


-2 


26 


+135 


59 


+0.61 







8354 


+118 


15 


+1 


81 


-128 


87 


-1 


66 


52743.0017 


0.2370 


-138 


96 


-0 


91 


+141 


72 


-2.22 


53158.8651 





2248 


-132 


66 


+4 


16 


+136 


23 


-6 


41 


52743.8681 


0.2599 


-138 


48 


-0 


26 


+136 


36 


-7.76 


53159.7571 





2779 


-137 


36 


-0 


97 


+141 


77 


-0 


38 


52745.9296 


0.6938 


+127 


96 


+1 


02 


-152 


37 


-13.76 


53182.6826 





3449 


-107 


94 


+7 


08 


+117 


19 


-2 


10 


52748.8870 


0.1854 


-124 


72 


+2 


87 


+133 


76 


+1.05 


53183.7214 





5713 


+55 


26 


-2 


71 


-66 


52 


+2 


40 


52751.9548 


0.8074 


+122 


44 


-4 


20 


-141 


35 


-3.07 


53183.8170 





6842 


+118 


57 


-5 


32 


-136 


77 


-1 


44 


52752.8104 


0.8176 


+118 


39 


-4 


95 


-137 


92 


-3.18 


53184.7284 





7602 


+139 


79 


+4 


76 


-147 


31 


-0 


01 


52769.7356 


0.8003 


+123 


65 


-5 


00 


-141 


11 


-0.66 


53189.7763 





7200 


+134 


61 


+1 


70 


-150 


90 


-5 


90 


52771.8062 


0.2449 


-139 


56 


-1 


14 


+137 


19 


-7.15 


53190.6836 





7912 


+129 


10 


-1 


72 


-146 


12 


-3 


34 


52773.8521 


0.6604 


+103 


84 


-10 


53 


-124 


10 


+1.06 


53192.6964 





1676 


-123 


51 


-2 


64 


+140 


40 


+14 


92 


52800.6851 


0.3407 


-116 


42 


+0 


57 


+128 


81 


+7.42 


53217.6764 





6602 


+113 


85 


-0 


42 


-123 


95 


+1 


11 


52802.6741 


0.6890 


+120 


03 


-5 


45 


-138 


94 


-1.90 


53452.8782 





3505 


-114 


77 


-2 


51 


+115 


01 


-1 


34 


52804.7941 


0.1920 


-127 


64 


+2 


03 


+145 


25 


+10.30 


53485.8782 





3118 


-121 


62 


+6 


74 


+125 


58 


-7 


96 



Note. — The O — C residuals correspond to the solution described in ^ 13.11 



TABLE 4 

Radial velocities for CV Boo from IPopperI (120001) . 



HJD 
-2 440 000 


Phase 


Star A 

km/s 


(O-C) 
A 


Star B 

km / s 


(O-C) 
B 


HJD 

-2 440 000 


Phase 


Star A 

km/s 


(O-C) 
A 


Star B 

km/s 


(O-C) 

B 


47198.0795 


0.6418 


+100 


59 


-4 


51 


-112 


85 


+2.43 


49117.9386 





3175 


-126.53 


-0.09 


+133.81 


+2 


32 


47254.9417 


0.7760 


+134 


29 


+0 


77 


-144 


43 


+1.25 


49202.6852 





3733 


-102.13 


-2.52 


+101.37 


-1 


49 


47397.6592 


0.2749 


-137 


69 


-0 


88 


+142 


48 


-0.12 


49204.7213 





7772 


+133.70 


+0.35 


-145.93 


-0 


43 


47695.7285 


0.1895 


-129 


31 


-0 


39 


+136 


69 


+2.54 


49204.7416 





8012 


+126.62 


-1.79 


-140.09 


+0 


10 


47696.7035 


0.3407 


-112 


58 


+4 


42 


+121 


75 


+0.35 


49204.7623 





8256 


+117.42 


-2.96 


-130.86 


+0 


69 


48080.7574 


0.7727 


+134 


39 


+0 


45 


-144 


60 


+1.53 


49496.9018 





7392 


+135.56 


+0.57 


-147.61 


-0 


36 


48081.6981 


0.8833 


+87 


15 


-3 


33 


-98 


55 


+0.90 


49496.9163 





7563 


+134.71 


-0.49 


-146.30 


+1 


18 


48312.0038 


0.7931 


+128 


21 


-2 


20 


-139 


74 


+2.60 


49583.6596 





1695 


-126.79 


-5.13 


+126.42 


+0 


08 


48344.9980 


0.7476 


+131 


13 


-4 


16 


-144 


15 


+3.42 


49583.6807 





1944 


-128.99 


+1.40 


+132.75 


-2 


98 


48345.8895 


0.8001 


+128 


12 


-0 


57 


-138 


79 


+1.70 


49907.6918 





7371 


+136.77 


+1.91 


-148.30 


-1 


19 


48345.9531 


0.8752 


+96 


80 


+1 


25 


-105 


62 


-0.74 


49907.7182 





7683 


+139.29 


+4.87 


-148.07 


-1 


43 


48819.6895 


0.1906 


-137 


81 


-8 


56 


+137 


28 


+2.78 


49907.7666 





8254 


+124.72 


+4.28 


-135.86 


-4 


24 


48819.7122 


0.2174 


-140 


25 


-4 


55 


+140 


79 


-0.64 


49907.7881 





8508 


+108.25 


-0.84 


-121.30 


-1 


88 


48819.7362 


0.2457 


-143 


18 


-4 


74 


+141 


52 


-2.85 


50176.9827 





6746 


+122.91 


+2.53 


-131.47 


+0 


11 


48820.6789 


0.3587 


-111 


27 


-3 


34 


+107 


70 


-4.03 


50177.0097 





7065 


+133.62 


+3.34 


-140.19 


+1 


99 


48820.7049 


0.3894 


-92 


10 


-2 


69 


+92 


66 


+0.66 


50177.9325 





7960 


+129.82 


+0.09 


-138.96 


+2 


65 


48822.7202 


0.7688 


+135 


49 


+1 


12 


-149 


18 


-2.59 


50177.9511 





8179 


+124.78 


+1.57 


-131.51 


+3 


09 


48822.8020 


0.8654 


+97 


14 


-4 


17 


-112 


48 


-1.42 


50177.9691 





8392 


+112.79 


-1.87 


-120.04 


+5 


37 


49116.9777 


0.1830 


-131 


43 


-4 


67 


+132 


23 


+0.41 


50177.9891 





8628 


+103.32 


+0.58 


-113.06 


-0 


46 


49117.8084 


0.1638 


-121 


42 


-2 


21 


+122 


85 


-0.85 


50178.0072 





8842 


+88.65 


-1.30 


-99.35 


-0 


46 


49117.8249 


0.1832 


-133 


33 


-6 


49 


+135 


33 


+3.42 


50178.0259 





9063 


+77.68 


+0.27 


-84.92 


-2 


06 


49117.8549 


0.2187 


-137 


56 


-1 


65 


+142 


79 


+1.14 


50616.7131 





8409 


+113.60 


-0.27 


-121.75 


+2 


80 


49117.8748 


0.2421 


-140 


97 


-2 


64 


+145 


24 


+0.99 



















Note. — Small corrections to the published values have been applied as described in ^ |2.3| The O — C residuals correspond to the 
solution described in ^ |3. 1| 
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TABLE 5 
Spot parameters for CV Boo. 







Co-latitude 


Longitude 


Radius 




Case 


Component 


(dcg) 


(deg) 


(deg) 


"^factor 


a 


pri 


8.062 


-5.47 


44.50 


0.82 






±59 


±90 


±50 


±5 


b 


sec 


7.926 


170.212 


47.76 


0.775 






±52 


±10 


±69 


±19 


c 


pri 


4.589 


13.058 


35.835 


0.5940 






±63 


±43 


±46 


±55 


c 


sec 


4.665 


151.43 


55 


0.881 






±43 


±38 


±1 


±11 


Note. 


— The spot 


co-latitudc is 


measured from the pole 


visible tc 



the observer, and the longitude is measured from the line joining the 
components' centers and increasing in the direction of orbital motion. 
The radius is measured as seen from the center of each component, and 
the temperature factor is relative to the unspotted photosphere. The 
uncertainties listed are in units of the last decimal place and correspond 
to the internal errors from the least-squares method. 



TABLE 6 

Light-curve solutions for CV Boo based on our V-band photometry. 



Parameter 



No spots Case (a) Case (b) 



Case (c) 



Parameter No spots Case (a) Case (b) Case (c) 



in 



Ob 

Teff ,B (K) . . . 

a. (R©) 

7 (km s~^ ) . . 
q = Mb/Ma- 



(^b/^a)v,0.25- 
cry (mag) 



87.651 
±42 

4.6752 
±76 

4.9014 
±78 

5632.8 
±1.6 

4.757 

±12 

-15.877 
±31 

0.9376 
±24 

0.741 
±10 
0.734 
0.0196 
0.428 
0.437 
0.715 
0.724 



86.891 
±34 

4.6591 
±48 

4.9614 
±50 

5628.1 
±1.4 

(4.748) 
(fixed) 

-15.889) 
(fixed) 

(0.9378) 
(fixed) 

0.746 
±4 

0.699 
0.0148 

0.429 
0.437 
0.715 
0.724 



86.650 
±33 

4.7844 
±49 

4.8186 
±47 

5656.1 
±2.4 

(4.748) 
(fixed) 

(-15.889) 
(fixed) 

(0.9378) 
(fixed) 

0.769 
±21 

0.829 
0.0147 
0.428 
0.435 
0.715 
0.722 



86.237 
±32 

4.7495 
±46 

4.8707 
±43 

5672.6 
±4.3 

(4.748) 
(fixed) 

(-15.889) 
(fixed) 

(0.9378) 
(fixed) 

0.778 
±14 

0.802 
0.0146 
0.428 
0.434 
0.715 
0.721 



^A,polc • ■ 
^A, point ■ 
'"A, side • ■ 
'"A, back • 
»'A,vol • • ■ 
'^B.polc ■ • 
''"B, point . 
''"B.sidc . • 
'"B.back ■ 
'"B.vol . • . 



PI- 



0.26533 
±67 

0.28105 
±90 

0.27032 
±73 

0.27726 
±82 

0.27189 
±74 

0.24054 
±98 

0.25177 
±117 

0.24410 
±104 

0.24935 
±111 

0.24483 
±105 
0.378 
0.390 



0.26646 
±49 

0.28268 
±67 

0.27167 

±53 

0.27877 
±62 

0.27252 

±55 

0.23699 
±86 

0.24757 
±101 

0.24042 
±91 

0.24535 
±94 

0.24107 
±91 
0.378 
0.390 



0.25798 
±47 

0.27188 
±61 

0.26254 
±51 

0.26870 
±58 

0.26327 
±51 

0.24573 
±84 

0.25825 
±101 

0.24972 
±90 

0.25548 
±96 

0.25049 
±90 
0.378 
0.388 



0.26028 
±46 

0.27478 
±62 

0.26501 
±50 

0.27142 

±57 

0.26577 
±51 

0.24247 
±81 

0.25423 
±97 

0.24624 
±86 

0.25168 
±93 

0.24697 
±87 
0.378 
0.386 



Note. — In all spotted solutions both components were assumed to rotate at a rate 1.0114 times faster than the orbital motion (see tj |3.1.1l l. 
For the no-spot solution the rotation is assumed to be synchronous. The linear limb-darkening coefficients [x) as well as the gravity-brightening 
exponents arc marked with an asterisk to indicate that they were changed dynamically du ring the iterations as T cff a nd log g changed. Th e 
gravity-brightening exponents varied over the mutually illuminated stellar surfaces following lAlcncar fc VazI 119971) and lAlencar et al 
and the values presented here are for the non-illuminated hemispheres. T^ff a was held fixed at 5760 K. The quantity (^b/^a) v,0.25 corresponds 
to the V-band light ratio at the first quadrature without considering the effect of spots, and (^b/^a)v is the mean light ratio outside of eclipse 
accounting for the spots and proximity effects. The uncertainties given on the left-hand side of the table (in units of the last decimal place) are 
the formal internal errors of the minimization procedure, while the ones on the right for the component radii account for the uncertainties of 
the gravitational pseudo-potentials as well as the mass ratio. The quantities rA,voi and rB,voi represent the "volume radius" for each star, i.e., 
the radius of a sphere with the same volume as the distorted stars. 



CV Boo 



TABLE 7 
Physical parameters of CVBoo. 



Parameter 



Absolute dimensions 

Mass (Mq) 

Radius (R0) 

log 9 (cgs) 

Measured t) sin j (kms~-^)... 

faaync sin j (kms'^) 

"sync sin j (kms"-"^) 

Radiative and other properties 

Teff (K) 

log L/Lq 

Mbol (mag) 

My (mag) 

Lb/La 

Distance (pc) 



Primary Secondary 



1.032 


± 


0.013 


0.968 


± 0.012 


1.262 


± 


0.023 


1.173 


± 0.023 


4.249 


± 


0.016 


4.285 


± 0.017 


73 




10 


67 


± 10 


78.5 




1.1 


72.7 


± 1.1 


77.6 


± 


1.1 


71.9 


± 1.1 


5760 


± 


150 


5670 


± 150 


0.197 


± 


0.048 


0.107 


± 0.049 


4.24 


± 


0.12 


4.46 


± 0.12 


4.32 


± 


0.12 


4.57 


± 0.13 






0.81 ± 0.13 








259 ± 16 





No te. — My an d Afboi were computed using bolometric corrections 
from iFloweil U996D along with Af®^^ ^ 4.732. The predicted asyn- 
chronous projected rotational velocities ^Jaaync sin i correspond to the 
values assuming the rotational period is Poj.b/1-0114 for both stars (see 
§[3T2j, while the Vsync sin z values give the result if P^-ot — ^orb- In both 
cases wc use the radius of the stars at quadrature. 



